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ABSTRACT 

Galaxy disks are shown to contain a significant population of atomic clouds 
of 100 pc linear size which are self-opaque in the 21 cm transition. These objects 
have H I column densities as high as lO^^cm"^ and contribute to a global opacity 
correction factor of 1.34±0.05 that applies to the integrated 21 cm emission to 
obtain a total H I mass estimate. High resolution, opacity-corrected images of the 
nearest external galaxies have been used to form a robust redshift zero distribu- 
tion function of H I, /{Nhi, X,z = 0), the probability of encountering a specific 
H I column density along random lines-of-sight per unit comoving distance. This 
is contrasted with previously published determinations of /{Nhj, X) at z = 1 and 
3. A systematic decline of moderate column density (18 < log(Nj^/) < 21) H I 
is observed with decreasing redshift that corresponds to a decline in surface area 
of such gas by a factor of five since z = 3. The number of equivalent Damped 
Lyman Alpha absorbers (log(N/^7-) > 20.3) has also declined systematically over 
this redshift interval by a similar amount, while the cosmological mass density in 
such systems has declined by only a factor of two to its current, opacity corrected 
value of (]gf^(z=0) = 5.4±0.9xl0-^ 

We utilise the tight, but strongly non-linear dependance of 21 cm absorption 
opacity on column density at z = to transform our high resolution H I im- 
ages into ones of 21 cm absorption opacity. These images are used to calculate 
distribution and pathlength functions of integrated 21 cm opacity. We suggest 
that that this z = calibration may also apply at higher redshift. In this case, 
the incidence of deep 21 cm absorption systems is predicted to show very little 
evolution with redshift, while that of faint absorbers should decline by a factor 
of five between z = 3 and the present. We explicitly consider the effects of H I 
absorption against background sources that are extended relative to the 100 pc 
intervening absorber size scale. Extended background sources result in dramat- 
ically altered distribution and pathlength functions which are insensitive to the 
predicted redshift evolution. Future surveys of 21 cm absorption will require very 
high angular resolution, of about 15 mas, for their unambiguous interpretation. 
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Subject headings: galaxies: individual (M31, M33, LMC) - galaxies: Local Group 
- galaxies: ISM - galaxies: evolution - cosmology: observations - radio lines: 
galaxies 



Introduction 



The 21 cm neutral hydrogen emission line has proven an extremely useful diagnostic 
of the interstell a r med ium of our own and nearb y galaxies since its first observation by 
Ewen &: Purcelll (1951) and MuUer &: Port ( 1951 ). And despite the fact that it was soon 
recognised ( iHagen &: McClairull954j ) that this transition could have significant opacity under 
interstellar conditions, it has often proven convenient to neglect that possibility when esti- 
mating the total H I mass of galaxies or the statistical properties of the neutral interstellar 
medium. Opacity estimation is straightforward in directions toward bright, compact back- 
ground continuum sources, b ut is considerably rn ore challenging to determine when based on 
the emission spectrum alone ( iRohlfs et al. 1119721 ) . A finite line opacity will cause systematic 
changes to the shape of emission spectra that can be used to estimate opacity for simple 
media and viewing geometries. Line-of-sight confusion encountered when viewing edge-on 
systems presents a very serious challenge to such an analysis given the likelihood for discrete 
features to overlap in both position and velocity. This has precluded detailed opacity mod- 
elling of the Galactic neutral hydrogen emission near and in the plane. Only by assuming a 
single universal effective ter nperature for al l the H I in the Ga l axy h ave some crude opacity 
corrections been made (e.g. ISchmidtl ( 1l957l ). [Henderson et al.l (119821 )). 



The first detaile d opacity modelling of the H I emission in an external galaxy has recently 
been undertaken by iBraun et al.l ( l2009l ) within M31. This was enabled by achieving 100 pc 
linear and 2 km s~^ spectral resolution at sub-Kelvin brightness sensitivity. High quality 
fits to individual spectra based on a physical model were possible over the majority of the 
disk, with the exception of those regions where the "warping" phenomenon of the outer disk 
leads to line-of-sight confusion. That study determined that opacity corrections could locally 
exceed an order of magnitude in the implied H I column density and that these resulted in a 
30% increase in the global H I mass estimate. But how representative are these values and 
to what extent are they influenced by the relatively high mean inclination, i = 78°, of the 
M31 disk or the finite spatial resolution? 



In this paper we extend the analysis of iBraun et al.l (120091 ) to the two other external 
galaxies for which the requisite physical resolution and sensitivity in the 21 cm H I emission 
line are available, M33 and the LMC. These galaxies have the virtue of a much more favorable 
disk orientation, with inclinations of about 50° and 25° respectively. The extreme proximity 
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of the LMC also permits linear scales as small as 15 pc to be probed with good sensitivity. 
Although small, this sample permits some general conclusions to be drawn regarding the 
statistical occurrence of opaque 21 cm H I emission and H I column densities. 

We begin with a brief description of the observations and their reduction in ^ and 
continue with their opacity analysis in ^ In ^ we present a statistical analysis of the 
occurrence of opacity corrected H I column density within all three galaxies and put these 
in a cosmological context by considering the prospects for detection of intervening 21 cm 
H I absorbers toward distant background source s. We assume a distance t o M31 of 785 kpc 



( iMcConnac 



lieetal 



2005 ). to M33 of 794 kpc flMcConnachie et al.ll2004f ) and to the LMC 
20091 ) throughout. We further a dopt a flat cosmolog ical model with 



of 51 kpc ( iKoerwer 
fl\ = 0.73 and a Hubble constant of 71 km s~^Mpc~^ (IHinshaw et al.ll2009[ ). 



Observations 



2.1. M31 



The observati ons and data reduction methods relevant to M31 are described in detail in 
Braun et al.l ( 120091 ). The most relevant aspects are the h igh quality tot a l pow er data obtained 
with the Green Bank Telescope (GBT) as described by lThilker et al.l (120041 ) combined with 
a 163 pointing mosaic obtained with the Westerbork Synthesis Radio Telescope (WSRT) in 
its "maxi-short" configuration. The data cube used for the spectral modeling analysis had 
angular resolution of 30" (114 pc), velocity resolution of 2.27 km s~^ and RMS noise level of 
1.0 K per channel. The relevant data attributes are summarized in Tabled] 



2.2. M33 



Total power observations of a 5° x 5° region centered on M 33 were carried out with the 
GBT during 2002 October using the same setup described by iThilker et al.l (120041 ) for the 
observations of M31. Interferometric data was acquired with the Very Large Array (VLA) 
in 1997, 1998 and 2001 under the proposal IDs AT206 and AT268. This consisted of a six 
pointing mosaic in the B and CS configurations and an additional 20 pointing mosaic in the 
D configuration. Some early results based on the B and CS configuration data were presented 
Thilker et al. (l2002h while an indep endent reduction of all three VLA configurations has 



m 



been presented by lGratier et al.l (120101 ). Joint deconvolution of all of the VLA pointing data 
was carried out after inclusion of the appropriately scaled GBT total power ima ges following 



the method employed for the reduction of the M31 data described in detail by iBraun et al. 
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( 120091 ) ■ The data cube used for the spectral modehng analysis had angular resolution of 



20" (77 pc), velocity resolution of 1.42 km s ^ and RMS noise level of 2.1 K per channel. 

2.3. LMC 

The observations and data reduction methods relevant to LMC are described in detail 



m 



Kim et al.l ( 120031 ). The most relevant aspects are the total power data obtained with the 
Parkes Telescope and its Multi-Beam 21 cm receiver combined with a 1344 pointing mosaic 
obtained with the Australia Telescope Compact Array (ATCA) in its 750A, 750C and 750D 
configurations. The data cube used for the spectral modeling analysis had angular resolution 
of 60" (15 pc), velocity resolution of 1.65 km s~^ and RMS noise level of 2.5 K per channel. 



2.4. Additional Targets 

While it would be very advantageous to extend our galaxy sample beyond the Local 
Group targets listed above, this is beyond the capabilities of current surveys and radio 
telescopes. 



For comparison, the THINGS survey (jWalter et al.l 120081 ) consists of H I imaging of 34 



relatively nearby galaxies using multiple configurations of the VLA. Although this survey 
achieves an angular resolution as high as 6 arcsec, which corresponds to less than 150 pc 
for the thirteen THINGS galaxies within about 4 Mpc, the corresponding RMS brightness 
sensitivity per velocity channel varies between about 10 and 20 K at this angular resolution 
while the velocity resolution is too coarse for several of the galaxies (5.2 km s~^) and only 
marginal (2.6 km s~^) for the remainder. As will become apparent below, parameterised 
spectral fitting is only found to be viable when at least one independent velocity channel 
exceeds lOo" and at least five exceed 5a. For this reason it is vital to achieve better than about 
2 K brightness sensitivity over 2 km s~^ at the required physical resolution of about 100 pc 
that matches intrinsic H I structure sizes. The THINGS survey has insufficient sensitivity 
by a factor of between 5 and 10 for this purpose. Remedying this shortcoming would require 
between 25 and 100 times longer integrations, or between 250 and 1000 hours per target 
with the EVLA. 

Another ongoing survey of ne arby galaxie s that might be considered of potential rele- 



vance is the HALOGAS survey of iHeald et al.l (1201 if ). That survey involves relatively deep 



(120 hours per target) observations of a galaxy sample within about 10 Mpc. Unfortunately, 
as that survey is designed primarily to detect diffuse gas in the extended environments of 
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galaxies it has only modest physical (about 1 kpc) and velocity (4.2 km s ^) resolution that 
make it unsuitable for parameterised spectral fitting of well-resolved structures. 



3. Spectral Analysis 



The method of spectral analysis has been described in detail in lBraun et al.l (120091 ). We 
repeat the most relevant aspects of that analysis below for clarity. 

Each line profile is modeled as a spatially resolved isothermal H I feature in the presence 
of turbulent broadening as 

TB{V)=Ts{l-exp[-T{V)]} (1) 

with 

,,,, 5.49 X 10-^^ Nhi f 

where the velocity dispersion, a, has units of km s~^ and is the quadratic sum of an assumed 
thermal and nonthermal contribution = {a'^ + aj^rp) with the thermal contribution given 
by ax = 0.093A/Tfc, for a kinetic temperature, T^. Such profiles have a Gaussian shape for 
low ratios of the H I column density (in units of cm~^) to spin temperature, Nhi/Ts, but 
become increasingly fiat-topped when this ratio becomes high. 

We found the best-fitting model spectrum for the brightest spectral feature along each 
line-of-sight when the observed peak brightness temperature exceeded 10 times the RMS 
noise and the truncated peak (where brightnesses greater than 50% of the peak were en- 
countered) spanned at least five independent velocity channels. The requirement for a well- 
sampled linewidth translates directly into one on the velocity resolution of the input galaxy 
data of better than about 2 km s^^. The truncation was done to both isolate single spectral 
components from possibly blended features as well as eliminating potential broad wings from 
the fit. In this way, the broad wings of secondary components along the same line of sight, 
for example from additional warm neutral gas with ~ 10^ K temperature, could be excluded 
from the fit. The data were compared to a pre-calculated set of model spectra spanning the 
range log(NH/) = 20.0 to 23.5 by 0.01, log(Tc;) = 1.2 to 3.2 by 0.01 and log{aNT) = 0.3 to 
1.5 by 0.04. A search in velocity offset was done over displacements of —4 to +4 by 1 km s~^ 
with respect to the line centroid as estimated by the first moment of each truncated peak. 

Only those fits with a normalized = ^^(A - Mif/{Na^) < 25 were retained; a 
cut-off based on the rejection level needed within M31 to exclude instances of line blending 
where the outer "warp" is seen in projection against the inner disk. In the case of M31 this 
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resulted in rejection of some 4% of the fits, while for both M33 and the LMC the fit rejection 
rate dropped to below 1%. 

Results of the spectral fitting were recorded as images of the physical parameters, 
N^f,Ts and cttvt together with the integral of brightness temperature over velocity that 
corresponds to the column density of the fit, J T^^^dV. This allowed calculation of a total 
corrected column density image from, 

AT™ = N^f + 1.823 X 10^*^ j Tl^'dV - j T^''dV^ (3) 

where J T'^°*dV is just the usual image of integrated emission. Those lines-of-sight with 
insufficient fit quality (x^ > 25) or insufficient peak brightness and width (at least one 
independent velocity channel exceeding lOcr and at least five exceeding 5cr) were assumed to 
have negligible opacity and simply assigned the optically thin value. 



4. Discussion 



4.1. Opacity-corrected images and the local 



We present images of the opacity corrected neutral hydrogen column density in Fig- 
ures [TH3J In contrast to images of the apparent column density, A''^^^ = 1.823 x 10^*^ / TsdV, 
which saturate below 10^^ cm~^ into an almost featureless fog, there is widespread fine struc- 
ture that extends to about 10^^ cm~^. It is interesting to note that comparably high H I 
columns have been observed in the Lyman Alpha absorption spectra of high redshift GRBs 
(IFynbo et al.ll2009l ). This contrast is illustrated in Figure H] where the corrected and apparent 
distributions are shown side by side for a region in the southeastern quadrant of the LMC 
that includes 30 Doradus (near (a, 5) = (05:40,— 69:00)). This figure also serves to demon- 
strate the exceptional detail provided by the LMC linear resolution of only 15 pc. While 
many of the high column density features in M31 and M33 are only marginally resolved 
with 100 pc resolution, the majority of the LMC features are well-resolved into structures 
that have coherent derived properties, despite the fact that every line-of-sight has been in- 
dependently fit with a model spectrum. The corrected total H I masses exceed those that 
follow from the assumption of negligible 21 cm opa city by facto r s of 1 .30, 1.36 and 1.33 for 
M31, M33 and the LMC respectively. As noted by iBraun et al.l (120091 ). the M31 value may 
be biased to a lower value since some portions of the disk are confused by the warp and 
could not be adequately modeled. Although our sample of galaxies with opacity-corrected 
H I column densities is very small, it does sa mple both the ma ss density peak and the 
highly populated tail of the H I mass function (jZwaan et al.ll2003l ) with apparent values of 
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Fig. 2. — Opacity corrected neutral hydrogen column density in M33. 
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Fig. 3. — Opacity corrected neutral hydrogen column density in the LMC 
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log{MHi) = 9.77, 9.13 and 8.61. And while all three galaxies are of moderat ely late type. 



the s ame is true of about 90% of galaxies within the complete HIPASS sample ( IZwaan et al. 



20031 ) both by their contribution to the total H I mass and in terms of space density. If 



these global correction factors are representative of galaxies in general, they suggest that the 
estimated redshift zero mass density, f^^j, determined from unbiased surveys of H I emission 
from galaxies in the local universe requires t his magnitude of u pward correction. Taking 
the av erage of the HIPAS S determined value (jZwaan et al.ll2003l ) and the recent ALFALFA 
value (iMartin et al.ll2010l ) after converting both to the assumed cosmology and scaling up 
by an opacity correction factor of 1.34±0.05 yields ^^^\{z = 0) = 5.6±0.9 x 10~^. 

The existence of additional gas, in exces s of that inferred from the published 21 cm and 
CO data has been previously suggested by [Bernard et al. J2008h in the case of the LMC 
from what has been termed "excess FIR" emission. Those authors have demonstrated that 
the dust optical depth is enhanced in regions of the highest apparent H I column density 
consistent with approximately twice the total gas mass being present. The opaque 21 cm 
emission that we document appears to account for some 35% of this FIR excess. This cool 
atomic hydrogen gas will likely be closely associated with diffuse molecular hydrogen gas 
that may be deficient in CO emission relative to the high density, self-gravit ating molecular 



clouds that are generally used to calibrate the CO to H2 conversion factor (e.g. iDickman et al. 



19861 ). A detailed analysis of atomic, molecular a nd dust tracers in th e LMC, including our 



opacity-corrected H I results will be presented in iThilker et al.l (j2012[ l. 



4.2. H I distribution and mass density functions 



It has been recognised for some time that images of 21 cm H I emission from nearby 
galaxies can be used to construct a distribution function that documents the probabihty to 
intercept a particular column d ensity on a random line-of-sight. The first application of this 
approach ( iRao fc Briggsl Il993[ ) employed the assumption that the apparent H I mass was 
distr ibuted smoothly over a ty p ical H I size that dep ends on galaxy type. Subsequent studies 



(e.g. iRyan- Weber et al.l (120031 ) . IZwaan et al.l (120051 )) have made use of interferometric images 
to estimate the areal coverage of specific column densities, but even these have been limited 
to a median physical resolution of some 1.4 kpc and no ability to correct for 21 cm opacity 
effects. Most recently, lErkal et al.l (120121 ) have used the THINGS galaxy images to this end, 
despite the fact that those data have only marginal resolution and inadequate sensitivity to 
enable resolved opacity corrections to be made, as discussed in §2.41 above. We have used 
the opacity corrected column density images (FigsilHS]) together with the s ensitive wide-field 
total power 21 cm H I emission data from the combined M31/M33 system ( iBraun fc Thilker 
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Table 1. Galaxy sample for spectral analysis 



Name 


Beam (arcsec) 


Beam (pc) 


ATb (K) 


(km s"^) 


Ref 


M31 


30 


114 


2.27 


1.0 


1 


M33 


20 


77 


1.42 


2.1 


2 


LMC 


60 


15 


1.65 


2.5 


3 



References. — (1) Braun et al. 2009; (2) Thilker et al. 2002; (3) Kim et 
al. 2003 




Right Ascension (JHOOO) Riglil Ascension (JHOOO) 

Fig. 4. — Opacity corrected (left) and apparent (right) neutral hydrogen column density in 
the southeastern quadrant of the LMC. 
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20041 ) and deep GBT imaging of the M31 environment (IThilker et al.ll2004j ) to construct the 
H I distribution function shown in Figure [5] that spans more than six orders of magnitude in 
H I column density from 10^^ - lO^^cm"^. 



The curves in Figure |5] were calculated from, 



f{Nni,X) 



Pa 



A(log(iVHIi)) 

/fo S,pr'" 'iVHiln(10)rflog(iVHi) 
where X represents unit co-moving depth along the line-of-sight, 

^0 



cm 



dX 



H{z) 



{l + zydz 



and 



Oi = ^,ln(10)(Mi/M,)"+^exp(-Mi/M,) 



(4) 



(5) 



(6) 



is the space density of galaxies with an apparent (i.e. uncorrected for opacity effects) H I 
mass Mi in units of Mpc^^dex^^ and A(log(A^Hii)) is the surface area subtended by H I in 
the column density interval (ilog(A^Hi) centered on A'hi for galaxy i. The area is measured 
by counting image pixels in each column density range and multiplying by the physical 
area of each pixel, pi = O jM] is the mass den sity due to galaxies of type i. The HIMF 



parameters were taken from lZwaan et al.l (120031 ): namely a faint-end slope, a = — 1.30±0.08, 
characteristic H I mass, log(M^,/M0) = 9.84 ± 0.06 and normalization, 6** = 7 .3 ± 1.8 x 10~^ 
Mpc~ ^dex~^, where Hq = 71 km~^Mpc~^ has been assumed throughout (IHinshaw et al. 



20091 ) and been used to rescale the HIMF parameters. The term, pav/^iPi in eqn. H] is 
introduced to normalize our measurements for three galaxies, to the global average density 
of apparentH I in galaxies, pav = 5.8 x ±1.0 x lO''M0Mpc~^, also taken from Zwaan et al. We 
stress that we make use of consistent quantities for this calculation, namely apparent galaxy 
H I masses and mass densities based on apparent mass. Alternatively, corrected quantities 
could have been used throughout. 

The solid curve in Fi gure El below log(Nfff) = 18.7 is taken from the wide-field (60° x 30°) 
WSRT total power data (IBraun fc Thilkerll2004[ ), the curve below log (N/:/f) = 19.8 from the 
extended (7° x 7°) GBT total power data of M31 (IThilker et al.ll2004l ). while the curve above 
log(Nj://) = 19.5 is that derived from the opacity corrected high resolution images presented in 
Figs. [IHSl For comparison, we also plot the distribution function without opacity correction 
as the dark dot-dash curve in the figure. The uncorrected data show an exponential upper 
cut-off that precludes column densities as large as lO^^cm"^. The opacity corrected data 
show an approximate power law extension to the distribution above about log(N///) = 21.5. 
We plot a power law with —3 index as the dot-dash line that provides a very good match to 
the measured distribution. Perhaps fortuitously, this power law index corresponds to what 
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Fig. 5. — Distribution function of neutral hydrogen column density (left) and the mass 
density function (right). Please see Figure [6] for greater detail of the inner portion of this 
figure. Solid curves represent our current, z = , determin a tion. Filled circles are used for 
the < z > =1 QSO absorption line dat a of iRao et al. I J2OO6I) and open circles for the 
< z > = 3 QSO absorption hne data of Prochaska fc Wolfel ( 2009 ) and O'Meara et al. 



(120071) . The dot-long-dash curves extend the < z > =3 data to lower Nuj (iProchaska et al. 
2OI0I ). The dot-dash curve is the z = determination of Zwaan et al. J2OO5I ) based on low 



resolution (1.4 kpc) H I emission images. The long-dashed and dashed curves represent 
the model we develop to represent the < z > =1 and 3 distributions. The short-dash, 
long-dash curve in the left panel is a —3 power law overlaid on the high Nm data. 
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would be expected for a population of resolved thin structures with some c haracterist i c face - 
on column density when viewed at random orientations, as first noted by iMilgroml ( 119881 ) . 
This may well be a plausible interpretation for observations such as presented here which 
resolve the atomic "skins" of what are likely to be randomly oriented clouds of molecular 
hydrogen on 100 pc scales within galaxy disks. 

Of particular note is the magnitude of the statistical errors in our determination of 
f{NHi,X). The solid curve above log(Nj|//) = 19.5 in the figure have been plotted three 
times: once as measured and again with both the addition and subtraction of the statistical 
la error calculated from the number of independent image pixels contributing to each column 
density bin. Only above log(Niif/) = 22.5 can these curves be distinguished from one another 
at all. The detailed shape of the combined curve above log(Nj://) = 19.5 is also reproduced 
in detail by the three galaxies individually, despite the fact that they represent different 
Hubble types and H I masses; ranging from SB(s)m, to SA(s)cd to SA(s)b and spanning 
an uncorrected log^Mni) = 8.6 - 9.8. From Figure 13 of IZwaan et al.l (120031 ) it is apparent 
that galaxies of these relatively late types contribute about 90% to the H I mass function in 
terms of both total mass and galaxy space density. The degree of inter-galaxy consistency 
is demonstrated with the error bars on the solid curve which represent RMS differences 
of the distribution functions of the three individual galaxies with respect to the combined 
fi^Hi, X). What hasn't been included in these error bars is the systematic error due to the 
uncertainty in the normalisation for galaxy space density, 6*^, (and hence pHi), of about 0.1 
dex. Errors in the normalisation would result in a vertical shift of the entire distribution 
function, but would not otherwise change the shape. 



We plot the parameterised fit to the H I distribution function of IZwaan et al.l (120051 ) 
as the short-dashed curve in Figs. [5] and [61 That curve is based on the analysis of a sam- 
ple of 355 galaxy images i n 21 cm H I emission made as part of the "WHISP" project 



(IVan der Hulst et al. Il200ll ). The distribution of Hubb le types in the WH ISP sample (IHolwerda et al, 
20 111 ) is similar to that seen in the complete HIPASS (IZwaan et al.ll2003l ) sample, consisting 
of about 90% galaxies later than Sb and 10% earlier. However, as noted at the outset of 
this section, it is limited by a median physical resolution of about 1.4 kpc and the inability 
to account for opacity effects. The authors suggest that their data might be representative 
in the range log(Nj|//) = 19.8 to 21, where sensitivity limits the low end and resolution lim- 
its the high end of that range. We plot the extr a polat ed extension of the fit beyond that 
range as a dotted line in the figure. IZwaan et al.l (120051 ) have explored the consequences of 
smoothing of the WHISP data to even coarser resolutions and note that this results in a 
steepening of the apparent fi^Nui^X) for log(Njy/) = 20 - 21. Conversely, they expect that 
achieving higher physical resolution would lead to a flattening of the derived f{NHi,X) in 
that column density range, which would then more closely resemble what we observe. It 
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should be noted that the ratio of median beam area between the WHISP sample and ours 
is a factor of about 400, implying that we sample about three times as many independent 
beam areas in our sample of only three galaxies as exist in the entire 355 galaxy WHISP 
sample. 

We have attempted to reproduce the effects of coarse physical resolution by smoothing 
the uncorrected (for opacity) images of our three galaxies to the median beam size of the 
WHISP sample, 1.4 kpc, and display these in Figure [71 As is clear from a comparsion of 
Figures [TH3] with Figure [TJ there is little useful information regarding the statistical ocurrence 
of intrinsic H I column densities at such a coarse physical resolution, given the extensive 
smearing of intrinsic peaks into larger areas of apparently lower column density. Although the 
combine d distribution function based on these three degraded images yields better agreement 
with the lZwaan et al.l (120051 ) curve, it still does not agree in detail. Looking at the fiN^i, X) 
contributions from the individual galaxies, it becomes apparent that the consequence of 
spatial smearing does not impact all galaxies to the same degree. While the redistribution of 
actual to apparent column densities has a moderate impact for the larger galactic systems, it 
becomes more dramatic for physically smaller ones. In the case of the LMC, the apparent area 
covered by log(Nj|//) = 20-21 gas is increased by about a factor of two by this smoothing. Its 
likely that this apparent enhancement of such column densities will be even more extreme 
in lower H I mass systems since they tend to have even smaller physical size. Since low 
H I mass systems strongly dominate total galaxy numbers, they will similarly influence the 
distribution function of a gala xy sample via their relative weighting, 6i in eqn. HI While we 
are not able to reproduce the IZwaan et al.l (120051 ) results in detail, it seems likely that the 
discrepancies are consistent with insufficient physical resolution in the Zwaan et al. study. 



Erkal et al.l (120121 ) have presented H I column density distribution functions based on 



the 34 THINGS galaxy images with physical resolution that varies between a bout 100 and 



500 pc and no ability to include opacity corrections. They compare these to the IZwaan et al. 



(120051 ) curve by arbitrarily normalizing the distribution from each THINGS galaxy to the 
Zwaan et al. curve at log(Nj^/) = 20, rather than utilizing the space density of each target 
galaxy. Based on this visual comparison they suggest that the H I distribution function is 
insensitive to physical resolution as coarse as several kpc. Our own comparison, discussed 
above and shown in Figs.Elandini utilizes an absolute normalization and clearly demonstrates 
that both high physical resolution and opacity corrections are vital in recovering the column 
density distribution function accurately from images of H I 21 cm emission. 

The filled circles with error bars in Figure El represent the < z > = 0.95 QSO 
absorption line data of iRao et al.l (120061 ). The open circles wi t h erro r bar s represent the 
high redshift QSO absorption line data of iProchaska fc Wolfe (l2009h and lO'Meara et al.l 



- 16 - 




Fig. 6. — Distribution function of neutral hydrogen column density (left) and the mass 
density function (right) as in Figure El but limited to a smaller range in H I column density. 




Right Ascension (J2000) „[gj,t Ascension (J2000) 



Fig. 7. — Apparent H I column density of M31 (left), M33 (center) and the LMC (right) as 
seen with an angular resolution tha t corresponds to 1.4 kpc, the median physical resolution 
of the "WHISP" survey used in the IZwaan et al.l (120051 ) study. 
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( 120071 ) ■ Those results are based on a sample of 738 damped Lya absorption systems (DLA, 
log{NHi) > 20.3) with < z > = 3.06 and 78 "Super" Lyman limit systems (SLLS, 19 < 
log(N//7) < 20.3) with < z > = 3.1. The < z > = 3.6 Lyman limit system data (LLS, 



17.5 < log(Nj://) < 19) from iProchaska et al.l (l2010l ) is plotted as the set of do t-long-dash 



lines. These parameterised fits, taken from Table 5 of IProchaska et al.l ( l2010l ). represent 
the range of allowed values. (We note in passing a minor typographical error in one of the 
tabulated parameters, the constant log k^^s ioi one of these curves should likely read —3.1 
and not —4.1.) While the QSO absorption line data is in good agreement with our z = 
distribution function for log(Nj|//) = 21 to 21.5, there is an interesting trend apparent for 18 
< log(Nj://) < 21, in which the areal coverage of such intermediate column density neutral 
gas declines continuously from z = 3 to z = 0. 

Another way of viewing all of the data presented above is to consider the corresponding 
mass density as function of H I column density. We plot this distribution in the right hand 
panel of Figure [51 The mass density was calculated from, 

Phi(A^hi) = mH—f{Nni, X)N^MW)dlog{Nni) (7) 

c 

where mn is the mass of an H I atom. This plot illustrates the range of column density 
that dominates the contribution to mass density. The trend noted previously for a declining 
contribution of modest column density gas (18 < log(N//7) < 21) with decreasing redshift 
is clear. That decline corresponds to about a factor of 5 since z = 3 at the DLA limit 
of log(Nj://) = 20.3, and about a factor of 2.5 relative to z = 1. The z = 3.6 LLS data 
demonstrate that this excess neutral gas disappears by log(Nj:/7-) = 17.5, where the areal 
coverage converges again with our z = distribution function. 

The plots of Figure [5] are presented again over a reduced range of column densities 
in Figure [6] to permit a more detailed comparison of the data and distributions, while the 
numerical values are listed below in Appendix A for convenience. 

Since the primary observed difference between the z = and higher redshift H I distri- 
bution functions is limited to an increase in the surface area of intermediate column density 
gas, we have considered a toy model for the z = 1 and z = 3 images that simply retains 
the smallest radii in the M31 image shown in Figure [H and replaces all structures at large 
projected radii (greater than 25 kpc or 1.8 degree), with the apparent column density version 
of the same image that is first made fainter by a factor of three (z = 1) and five (z = 3) and 
scaled up in linear size by a factor of 1.55 (z = 1) and 2.3 (z = 3). This could represent 
an extended intermediate column density halo around the high column density disk, that is 
larger in area than the actual z = galaxy by a factor of 2.4 for z = 1 and 5.1 for z = 3. The 
resulting distribution function and mass density function are overlaid on Figs. and E] as 
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the long-dashed and short-dashed curves. These curv es provide quite a good representation 
of the QSO absorption hne data at z = 1 and z = 3 of iRao et al.l (l2006l ). iProchaska &: Wolfe 
fl2009h and lO'Meara et al.l fcoOTh . 



While this toy model is not physically motivated, it serves to demonstrate the types of 
structures needed to account for the enhanced contributions of intermediate column density 
gas in the high redshift data, and achieves these with only minimal changes to the z = 
images. Such a distribution might arise if gas accretion rates fror n the intergalactic medium 
were substantially enhanced in the past. The recent simulations of lKeres et al.l ( 120091 ) suggest 
that the total gas accretion rate was enhanced by factors of about 3 and 10 at redshifts of 1 
and 3 over that at z = 0, with the enhancement of "cold mode" accretion of mildly ionized, 
10^~^ K gas, being even more dramatic at these redshifts. 

An alternative model, that is perhaps more plausible, might consist of a population of 
dwarf satellites surrounding each massive high redshift galaxy, the majority of which will 
merge with the parent by z = 0. Since the spatial concentration of such a system will 
be similar to the halo model and the areal coverage of both high and intermediate column 
density gas will need to be identical (to match the QSO absorption line data) there will likely 
be only minor practical differences between these two scenarios for the purpose of a statistical 
description. We have not attempted to pursue modeling of a merging satellite population at 
this time because of the many additional free parameters that would be introduced. 



4.3. H I pathlength functions 



A quantity that is often more useful in planning and interpreting surveys is the path- 
length function, dN/dX, the number of intervening systems per unit co-moving distance 
along any line-of-sight that exceed a particular column density. This is formed by integrat- 
ing the positive tail of the H I distribution function and is shown in Figure El The two 
model curves discussed above to represent the z = 1 and 3 QSO absorption line data are 
also plotted in the figure. The total number of equivalent Damped Lyman Alpha (DLA) 
absorption systems, defined to be those with H I column density, log^Nui) > 20.3, is 
nDLA{X, z = 0) = 0.026±0.003 for our redshift zero deterrn ination. This is almost a factor of 
two lower than the value estimated by IZwaan et al.l (120051 ). nDLA{^,z = 0) = 0.045±0.006, 



and seems a consequence of the insufficient physical resolution (a median of 1.4 kpc relative 
to a structure scale of 100 pc) of that stud y. Our me a suren aent i s compared with the values 
of riDLAiX, z) published in the studies of iRao et al.l ( 12006! ) and jProchaska fc Wolfd (120091 ) 
in the center panel of Figure m There has been a continuous decline in nDLA{X) with red- 
shift as noted by many previous authors. Our redshift zero measurement demonstrates an 
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ongoing decl i ne sin ce z = 1, rather than flattening to a constant value as suggested by the 
Zwaan et al.l ( 120051 ) point. We include the unLAi^, z) values of our model for the z = 1 
{riDLAi^, z = 1) = 0.042) and z = 3 {unLAi^, z = 3) = 0.064) distributions on the plot 
to demonstrate the degree to which they are in agreement with the published data. The 
right-hand panel of Figure [8] presents the redshift evolution of the mass density in H I 

due to systems with column densities Ioq^Nhi) > 20.3. This was calculated for redshift zero 
by considering the integral mass above this threshold, which amounts to 96% of the total 
H I mas s in galaxies or Q ^^^jz = 0) = 5.4±0.9 x 10~^. The data poi nts in the figure are 
those of iRao et al.l (120061 ) (after rescaling to an H I only value) and of iProchaska fc Wolfe 
(120091 ). For comparison we also plot the mass densities in our redshifted model distributions 
to demonstrate their consistency. 



4.4. Relating 21 cm absorption to H I column density 



The 21 cm H I transition offers great promise for documenting the evolution of neutral 
gas content with cosmic time. However, there are specific issues which need to be overcome 
before it can fulfill that promise. In the case of 21 cm emission, the issue is sensitivity. Even 
the deepest integrations with existing fac ilities have only been a ble to push back this frontier 
to z ~ 0.2 for i ndividual target detection (jVerheijen et al. Il2010l ). And while so-called "stack- 
ing" analysis (ILah et al.l 120071 ) can provide some statistical information on gas content at 
larger distances, the interpretation of such results can not be disentangled from the definition 
of the sample of objects to be co-added. Reaching down to below M^, on individual targets 
at redshifts exceeding unity will require about 100 times the sensitivity that is currently 
available, which can only be achieved with a hundred-fold increas e in collecting area such a s 
planned for Phase 2 of the international Square Kilometre Array (jCarilli fc Rawlingsi 120041 ). 
For 21 cm line absorption, the major limitations are firstly, the need for a large instanta- 
neous field-of-view which would enable effective blind surveys toward background continuum 
sources and secondly, an effective calibration of integrated absorber strength with the asso- 
ciated H I column density. The limitation of field-of-view is currently b eing addressed with 
Phased Array Feed (PAF) technology deployed on parabolic reflectors (IDeBoer et al.l 120091 ) 
that will provide as much as 30 deg^ of instantaneous sky coverage. The limitation of 21 cm 
absorber calibration is one that we will attempt to address here. 

The necessity for calibration of the 21 cm absorption line strength is determined by 
the intrinsic temperature sensitivity of the 21 cm H I line opacity, which varies inversely 
with temperature, as already noted in eqn. [2] as well as on possible nonthermal contribu- 
tions to the linewidth. Since any particular line-of-sight is likely to traverse more than one 
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temperature of medium, there might be no simple correspondence of the 21 cm hne opacity 
with the underlying H I column density. Fortunately, there appear s to be a parti c ularly 
good one-to-one correspondence of these quantities observationally. iKanekar et al.l (1201 if ) 
have recently demonstrated the remarkably tight relationship between integrated 21 cm 
absorption opacity and the apparent (ie. uncorrected for opacity) column density deter- 
mined from complimentary 21 cm emission observations of neighbouring lines-of-sight. We 



present a simi 



Kanekar et al- 



ar plot in Figure [9] where the Galactic 21 cm absorption and emission data of 
(120111) are plot t ed as filled circles, together with similar d ata obtained toward 



M31 b y iBraun &: WalterbosI (Il992l ) as open circles and the LMC by iMarx-Zimmer et al. 



( I2OOOI ) as open squares. Only the high signal-to-noise absorption detections (exceeding 4cr) 
are plotted. The curves in the figure represents the expectation for a "sandwich" geometry 
of cool gas wit h properties (Tr , jr) su rrounded by layers of warm gas with properties (T^, r^) 
as discussed in iKanekar et al.l (l201l[ ) from. 



N'hi 



Noe-< + N^{l-e-< 



for a measured "apparent" column density, A^^j (based on the integral of observed brightness 
temperature), a threshold column density (where Tc 0), A^'o ~ T^r^„Al^, a saturation 
column density (where Tc 00), A^oo ~ (^c + TuiT^/^^/SV and an effective opacity, r^. The 
effective opacity is related to the measured integrated opacity by the effective linewidth, A\^, 
as / rdV = r^AV". The solid curve is for (Nq = 1.25 x 10^° cm-^, = 7.5 x 10^^ cm'^, 



/SV = 15 km/s) while the dashed and dot-dash curves are for (No 



l20 



5.0x10^1 cm-2, AV 
respectively. 



20 km/s) and (Nq = 2 x 10^° cm'^, Noo = 10^^ cm'^, AV = 10 km/s) 



The "sandwich" geometry is meant to represent the role of a "temperature-shielding" 
column density of the Warm Neutral Medium (WNM) that absorbs local EUV and X-ray 
radi ation and is a requir ement for condensation of the Cool Neutral Medium (CNM) within 



Wolfire et al.l 120031 ). It should be noted that this "temperature-shielding" column is 



(e.g. 

quite distinct from the "self-shielding" column density of only a few times 10^^ c m~^ that 



mark s the onset of recombination of hydrogen within an ionizing radiation field (e.g. iDove fc ShuU 



19931 ). Once neutral columns of about 10 cm are achieved, the expectation is that the 



ionization fraction will be minimal, even in extreme radiation fields. Such a model provides 
a very good characterisation of the observable properties of the absorbing medium in the 
Galaxy, where it is best constrained by data, but also of the available extragalactic data 
from M31 and the LMC, in which both higher and lower metallicity and radiation field en- 
vironments are probed. The dashed and dot-dash curves provide a reasonable estimate of 
the scatter observed in this relationship. Since the curves drawn in Figure [9] are based on a 
physical model, we can also write the corresponding equation for the actual rather than the 
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apparent H I column density, 

Nhi = No + {N^ - -fK (9) 

The curves based on the same parameters as previously are drawn in the right-hand panel 
of Figure [9l While these are unchanged at low column density, they do not exhibit the 
saturation effect at high columns that afflict observations of 21 cm H I emission. 



4.5. 21 cm opacity distribution function 

With the column density calibration of 21 cm absorption line strength of the previous 
sub-section in hand it becomes possible to construct the distribution function of 21 cm 
opacity. Eqn. [Hlwas used to transform the three images of opacity-corrected column density 
(Figs. [IH3]) to corresponding images of 21 cm absorption line strength. The transformation 
was done using the nominal transformation depicted in Figure [9] with parameters (Nq = 
1.25 X 10^° cm~^, Noo = 7.5 x 10^"*^ cm~^, AV = 15 km/s) as well as the two envelope curves 
with (No = 10^0 cm-2, N^o = 5.0 x lO^^ cm^^ = 20 km/s) and (Nq = 2 x 10^° cm-^, 
Noo = 10^^ cm~^, AV = 10 km/s) in order to track the intrinsic scatter and uncertainty. The 
distribution function was then calculated from eqn. IHbut with J rdV in the place of Nhi. 
This is shown in Figure [10] with the solid curve illustrating the nominal z = relation, while 
the dashed and dot-dashed curves span the scatter and uncertainty. What is notable is the 
flatness of the distribution below about J rdV = 0.1 km s^^. The "temperature-shielding" 
threshold implies that there will be essentially no additional systems seen with an intrinsic 
absorption depth below this value. 

Distribution functions of the 21 cm absorption for the toy models developed above to 
describe the z = 1 and 3 QSO data are also plotted in Figure [TD] as the dotted and short- 
dashed curves utilizing the nominal z = transformation. There is considerable uncertainty 
in the applicability of the z = relation at these higher redshifts. The intergalactic ionizing 
radiation field is expected to increase dramatically with redshift in concert with the evolution 
of the massive star formation rate and active galaxy incidence. While this will have a major 
influence on the "self-shielding" column density (of a few times 10^^ cm~^) for the onset of 
the ionized to neutral transition at the galaxy-IGM interface, it is less clear what the impact 
might be on the "temperature-shielding" column (of ~ 10^° cm~^) for the warm to cool 
transition that occurs deep within individual galaxy disks. This is b ecause the radiatio n 



field in this latter case is strongly dominated by local contributions (e.g. IWolfire et al.ll2003l ). 
It is likely that the intrinsic scatter of the z = relation in Figure [H] is already a consequence 
of the large variations in local radiation field intensity that naturally occur within individual 
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Fig. 8. — Pathlength function of H I column density (left), as well as redshift evolution of 
the total number (center) and H I mass density (right) of DLA systems. 




Fig. 9. — The relation of integrated 21 cm absorption with apparent neutral hydrogen column 
density (left) and total neutral hydrogen column density (right). 
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galaxies. Future work will be vital to test this hypothesis in more extreme environments, 
such as the disks of nearby star-bursting galaxies. We will continue with our discussion of 
possible evolutionary trends of absorption properties with this caveat in mind. 

Since deep 21 cm absorption features, / rdV > 5 km s~^, require high column density, 
Nhi > 10^^'^ cm~^, lines-of-sight and the incidence of such column densities does not appear 
to evolve significantly with redshift (Figure [6]), there is also no evolution in the distribution 
function of 21 cm absorption in Figure [TD] for such deep features. Redshift evolution is 
confined to the relatively weak 21 cm absorbers with strength J rdV < 1 km s~^ that 
correspond to Nhi < 10^^ cm~^. 



4.6. The pathlength function of 21 cm absorption 

Just as was done above for the H I column density, it is possible to calculate the 
corresponding pathlength function, dN/dX, the number of 21 cm absorption systems per 
unit distance that exceed a particular absorption depth, J rdV. This is formed by integrating 
the positive tail of the 21 cm absorption distribution function and is shown in Figure [TT] with 
the same curves used previously to depict the z = 0, 1 and 3 predictions. The same trends 
noted previously for the distribution function apply here as well, namely an absence of 
redshift evolution for deep absorbers and number evolution by a factor of several for the 
faint absorbers. 



4.7. The impact of finite background source size 

The preceeding discussion has assumed that the background source used to probe inter- 
vening absorbers has a small angular size relative to the most compact structures within the 
absorbing system. What if this is not the case? Since we already have images that repre- 
sent absorption depth, it is straightforward to simulate the observable effect of an arbitrary 
background source size and structure. This can be done by evaluating the convolution of the 
absorption depth image with an image of the background source, normalized to preserve the 
integrated absorption depth. The distribution function of such convolved images can then 
be calculated in the same way. This has been done in the right hand panels of Figs. [10] and 
[TT] for an assumed Gaussian background source of 9 kpc FWHM. This choice of background 
source model is motivated by the 1.3 arcsec median source size observed for mJy brightness 



radio sources at GHz frequencies (IWindhorst II2003I ) which subtend more than 9 kpc for all 



redshifts greater than about 0.7. This should be regarded as a "worst case scenario" from 
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Fig. 10. — Distribution function of integrated 21 cm absorption at full resolution (left) and 
observed with 9 kpc effective resolution (right). 




Fig. 11. — Pathlength function of integrated 21 cm absorption observed at full resolution 
(left) and observed with 9 kpc effective resolution (right). 
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the standpoint of absorption signal dilution, since there is some probability that a compact 
component may be present, even in an intrinsically extended background source. 

The impact of a significantly extended background source on the 21 cm absorber dis- 
tribution and pathlength functions is dramatic. As can be seen in the right hand panels of 
Figs. [To] and [TTl deep absorption systems will not be observed at all in this instance, while 
an enhanced number of apparently faint absorbers are seen instead. This redistribution of 
intrinsically deep into apparently faint absorbers also has as a consequence that there is no 
predicted observable evolution at all in the 21 cm absorption properties with redshift. 

As noted above, this assessment must be regarded as a "worst case scenario" since the 
median source size has been distributed smoothly in our simulation. Better insights into 
the typi cal structures seen in rad io sources at GHz frequencies are provided by the recent 



study of iMiddelberg et al.l ( 120111 ) who studied a sample of nearly 100 moderate brightness 
radio sources with angular resolution of about 20 mas. Overall, some 20% of these sources 
were found to have a substantial (50 - 100%) flux contribution from unresolved structures. 
However, when considering sources fainter than 1 mJy at 1.4 GHz, this fraction declined to 
only 7%. 

What the dramatic contrast between the left- and right-hand panels of Figs. [10] and [11] 
underlines is the importance of isolating the flux contribution of suitably compact structures, 
subtending less than about 100 pc at the distance of the absorption system, in interpreting 
any measurement of 21 cm absorption. Interpretation is essentially impossible without this 
information. The implication for large surveys of H I 21 cm absorption at redshifts z > 0.5 
is that high resolution imaging, with about 15 mas beam size, will be needed for the useful 
interpretation of any detections, either during the survey itself or as follow-up. At the 
relevant observing frequencies, v < 700 MHz, this implies a need for baselines greater than 
about 5500 km. 



5. Conclusions 

The analysis of high resolution and sensitivity H I datacubes of Local Group galaxies 
has permitted some fundamental conclusions to be drawn: 

1. Galaxy disks contain a significant population of H I features which are self-opaque in 
the 21 cm transition. Such features are intrinsically compact, with narrow dimensions 
that are typically only about 100 pc, which determines the resolution required for their 
study. Such features can locally disguise the presence of high atomic column densities 
in the range 22 < log{NHi) < 23, that are not apparent in images of the integrated 
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21 cm linestrength. Globally, such features account for about 34±5% more atomic gas 
than is derived under the assumption of negligible self-opacity. 

2. Application of this opacity correction to recent estimates of the cosmological mass 
density of H I within local galaxies yields, ^h\{z = 0) = 5.6±0.9 x 10~^ assuming that 
our limited galaxy sample provides a representative determination. 

3. The resolved, opacity corrected images of Local Group galaxies permit calculation of 
a robust redshift zero H I distribution function, f{NHi,X), with excellent statistical 
properties (based on lOOO's of independent sight-lines) and a high degree of internal 
consistency amongst galaxy types between Sb and Sm. While similar to high redshift 
determinations, there is evidence for a systematic decline in intermediate column den- 
sity gas, 18 < log{NHi) < 21, that can be well modelled by a decrease in the surface 
area of galaxy halos (or associated satellites) by a factor of about 2.5 since z = 1 and 
a factor of about 5 since z = 3. 

4. The number of equivalent Damped Lyman Alpha (DLA) H I absorbers per unit comov- 

ing length, n£)/^^(X, ^ = 0) = 0.026±0.003, is significantly lower than seen at higher 
redshifts by factors as large as about four and follows a smoothly declining trend. Sim- 
ilarly, the cosmological mass density in DLAs, Ogf^(^ = 0) = 5.4±0.9 x 10-^ follows 
a smoothly declining trend, albeit one that sees only a factor of two decline between 
redshift 4 and the present. 

5. There is a tight, non-linear correlation between the observed 21 cm absorption opac- 
ity and the integrated 21 cm emission brightness from adjacent sight-lines through 
nearby galaxy disks which permits "calibration" of our H I column density images into 
corresponding ones of integrated 21 cm opacity. 

6. We present predictions for the distribution and pathlength functions of 21 cm inte- 
grated opacity at both redshift zero and at high redshift for the case of no evolu- 
tion in the local opacity calibration. The prediction is that faint 21 cm absorbers, 
J rdV < 1 km s~^, have declined by a factor of about 2.5 since z = 1 and a factor of 
about 5 since z = 3, while the number of deep absorbers, J rdV > 5 km s~^, has not 
evolved. 

7. We explore the impact of a finite background source size on the distribution and path- 
length functions of 21 cm integrated opacity. Adoption of the observed 1.3 arcsec 
median source size of mJy radio sources dramatically changes the apparent opacity 
distributions and effectively precludes any physical interpretation of the measurements. 
Future surveys of 21 cm absorption will require a 15 mas beamsize for their unambigu- 
ous interpretation. 
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A. Numerical values of the redshift zero H I distribution function 

The numerical values of the H I distribution function presented in this paper arc tab- 
ulated below. The data source for the different column density ranges are indicated. The 
error estimates in the third column arc determined from the number of indcpcndant image 
pixels used in the evaluation for the data of references 1 and 2, while for reference 3 this is 
determined from the RMS deviation amongst the three contributing galaxies, which exceeds 
the statistical error estimate by factors of between four and 20. 
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Table 2. HI distribution function at z = 



7 /AT" \ r —91 

log{NHi)[cm ^] 


log{f{NHi,X) 


^log{f{NHi, X) 


Ref 


16.75 


-17.28 


0.093 


1 


17.00 


-17.43 


0.083 


1 


17.25 


-17.74 


0.089 


1 


17.50 


-18.38 


0.140 


1 


17.75 


-18.82 


0.174 


1 


18.00 


-19.25 


0.213 


1 


18.25 


-19.69 


0.265 


1 


18.50 


-20.14 


0.334 


1 


18.75 


-20.44 


0.353 


1 


18.58 


-20.30 


0.188 


2 


18.84 


-20.60 


0.195 


2 


19.11 


-21.06 


0.246 


2 


19.37 


-21.40 


0.269 


2 


19.63 


-21.82 


0.322 


2 


19.89 


-22.15 


0.348 


2 


19 62 


-21.91 


321 


(J 


19 80 


-22.10 


227 




1 9 97 


-22 29 


121 




20 15 


-22.43 


027 


(J 


20.32 


-22.53 


0.032 


3 


20 50 


-22.64 


058 


3 


on {^'7 


OO 70 

-zz. / y 


v.OOi 


Q 
O 


20.85 


-22.95 


0.044 


3 


21.02 


-23.05 


0.109 


3 


21.20 


-23.19 


0.182 


3 


21.37 


-23.47 


0.086 


3 


21.55 


-23.92 


0.084 


3 


21.73 


-24.47 


0.139 


3 


21.90 


-25.16 


0.161 


3 


22.07 


-25.76 


0.195 


3 


22.25 


-26.30 


0.152 
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Table 2 — Continued 



log(Nni)\cm-^] 


log{f(Nni.X) 


Mog(f(N,ii.X) 


Ref 


22.43 


-26.73 


0.161 


3 


22.60 


-27.09 


0.315 


3 


22.77 


-27.52 


0.332 


3 


22.95 


-28.25 


0.461 


3 
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